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We present a general theory of electric field effects in liquid crystals where the dielectric 
tensor depends on the orientation order. As applications, we examine (i) the director fluc- 
tuations in nematic states in electric field for arbitrary strength of the dielectric anisotropy 
and (ii) deformation of the nematic order around a charged particle. Some predictions are 
made for these effects. 
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1. Introduction 

Phase transitions occur in various systems in electric field. To treat such problems we 
need to construct a Ginzburg-Landau free energy including electrostatic interactions and 
electric field effects under each given boundary condition.^) The gross variables include the 
order parameter ip, the polarization p, and the charge density p. (Note ijj = p m. ferroelectric 
systems.) In this letter we present such a theory for liquid crystals in nematic states, since 
there seems to be no systematic theory of electric field effects for anisotropic fluids. We also 
examine the deformation of the director around a charged particle, since the nematic order 
around a neutral colloid particle or emulsion droplet has been extensively studied.^' We 
shall see that the charge-induced orientation is intensified with decreasing the particle radius 
R and/or increasing the charge number Z, while the surface anchoring of a neutral particle 
can be achieved for large radius because the penalty of the Frank free energy needs to be 
small. 

2. Liquid crystals between capacitor plates 

In liquid crystals the order parameter ij: represents the orientation tensor near the 
isotropic-nematic transition or the director n in nematic states. We also use the vector nota- 
tion p to represent the set {pa\ for the ion densities pa (a = 1, 2, • • • ). We divide the total 
free energy functional F = F^h + Fgt into a chemical part Fch and an electrostatic part Fgt- 
Here, 

The first contribution Fq is assumed to be independent of p, but there arises a coupling 
between ij) and p in the presence of flexoelectricity.^'^-' The tensor x*"' is the inverse matrix 
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of the the electric susceptibility tensor Xij ■ We shall see that the local dielectric tensor eij is 
related to Xij ^ 

Eij = 6ij +4Trxij- (2) 

As is well-known, Sij depends on the orientation order. In particular, in nematic states the 
following form has been assumed in agreement with experiments:^^ 

Eij = EoSij + EiniUj. (3) 

The coefficients eq and Ei satisfy eo > 1 and eq + £i > 1 because Xij should be a positive- 
definite matrix from the thermodynamic stability. 

Next we consider the electrostatic part Fgt- A typical experimental geometry is shown 
in Fig. la, where wc insert our system between two parallel metallic plates with area S and 
separation distance L. We assume S^/^ » L and neglect the effects of edge fields. The z axis is 
taken perpendicularly to the plates. Let the average surface charge density of the upper plate 
be (Tex and that of the lower plate be — Cex- The total charge on the upper plate is Q = Saex- 
The electric potential (j) satisfies = at the bottom z = and = $ at the top z = L, 
where $ is the potential difference between the two capacitor plates. The electric induction 
D = E + Anp satisfies 

V D = -V'^4> + 47rV • p = Attp, (4) 

where p = eZ^Pa is the charge density with Z^e being the charge of the ion species a. 
The boundary conditions at z = and L are Ex = Ey = and Dz = —Anae^.^^ With these 
relations Fgi is of the form, 

F,t = J drE^Sir. (5) 

In fact, if infinitesimal deviations dp, dp, and SQ are superimposed on p, p, and Q, the 
incremental change of Fgt is given by 

dF,t = J dr[-E ■ 5p + (l)5p\ + ^5Q, (6) 

where use is made of 5{E'^) = -2V0 ■ 5D - 8nE ■ 6p = -2V{<j)5D) + 8Tr{(t)dp -E-Sp). The 
right hand side of (6) represents the increase of the electrostatic free energy due to the small 
variations of p, p, and Q. 

If we minimize F with respect to p for fixed tp, p, and Q, we require SF/Sp = 0. Then 
(1) and (6) give 

p =x -E, D =e -E, (7) 

which indeed yield (2). Then we may rewrite F as F = Fq + F^ with Fe = J dr X^^PiPjl'^'^ 
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(8) 



Under these relations and using 5x^^ = ~J2ke X^'^X^^^XU we obtain 



dr 



5F = 6Fo + 

Prom (4) and (7) the electric potential satisfies 

-V- s -Vcf) 



(t)6p - ^E -Sl -E 
Stt 



(9) 



(10) 



Thus 4> and E = — V0 depend on the orientation order via (3) and on p under each boundary 
condition. In the literature,^) however, the electric field E in (9) has been replaced by the 
average constant field Eq, where F = Fq — £i J dr{Eo ■ n)^/87r +const.+0(ef ). This is the 
first order approximation with respect to ei valid for |ei| <^ eq. If ei + £o is considerably 
smaller or larger than there can be unexplored regimes of strong polarization anisotropy. 

So far the capacitor charge Q is a control parameter and the potential difference $ is a 
fluctuating quantity dependent on p and p. We may also control ^ by using (i) a battery at a 
fixed potential difference or (ii) connecting another large capacitor in parallel to the capacitor 
containing our system as in Fig. lb. As will be discussed in the appendix, the appropriate free 
energy functional is given by the Legendre transformation,^^ 



G = F-^Q. 

Prom (9) the incremental change of G is obtained as 



6G = SFo + J dr 
3. Director fluctuations in nematic states 



(bSp - ^E -Se -E 
Stt 



(11) 



(12) 



As a first application of our theory, we apply an electric field Eq to a nematic liquid 
crystal considerably below the transition without charges {p = 0), where the orientation 
order is represented by the director vector n. Here we do not perform the expansion with 
respect to £i. For simplicity, we assume that n is along the z axis on the average under 
the homeotropic boundary condition and Si > 0. If Si < 0, the orientation along the z axis 
becomes unstable with increasing Eq above a critical value (oc L~^).^) Then n is written as 
Bz + Sn with 6nz = — [{Snx)"^ + (5%)^]/2, where is the unit vector along the z axis and 6n 
is the deviation nearly perpendicular to the z axis. By solving (10) under (3), we may expand 
the electrostatic free energy Fe = F — Fq in powers of Sn. The fluctuation contributions on 
the bilinear order are written as 



(13) 
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where q = q^^q represents the direction of q and n{q) is the Fourier transform of 6n. The 
wave number |g| is assumed to be much larger than the inverse system width L~^. The large 
scale fluctuations are omitted in (13). 

Expressing Fq in terms of the Frank constants Ki, K2, and K^, we obtain the director 
correlation functions, 

' (<5n,(g)^n,(qr) = ^ + ^fl-lV (14) 



where i,j = x, y, qj_ = ql + q^, and 

n = K^ql + Kiql + ei{eQ + £i)(£o + eiqly^E^/An, (15) 

r2 = K^ql + K2ql + £iEI/Att. (16) 

If £1 > 0, the correlation length ^ is given by ^ = {Ai^K / £1)^/"^ E^^ , where K represents the 
magnitude of the Prank constants The coefficient ri depends on q even in the limit (7^0, 
which is not the case in the previous literature. The scattered light intensity is proportional 
to the following,^) 



kBTel r2 \ri r2 J q\ 

where i and / represent the initial and final polarizations. The vector a is defined by 

In addition we examined anisotropy of the turbidity (form dichroism) in a nematic state under 
electric field in Ref.[l]. 

4. Orientation around a charged particle 



We place a charged particle with radius R and charge Ze in a nematic state, where n 
is aligned along the z axis or n ^ Bz fax from the particle. Let the density of such charged 
particles be very low and its Coulomb potential be not screened over a long distance A. Prom 
(3) and (9) the free energy change due to the orientation change Sn is given by 



6F = - dr 



KV^n + -^(E-n)E 
47r 



• 6n, (19) 



where we have assumed the single Prank constant K {Ki = K2 = = if), so Fq = 
/ drK\\/n^ /2. If the coefficient £1 is considerably smaller than £0, the electric field E near 
the particle is of the form — (Ze/£or^)f and the electrostatic energy F^ is approximated by 



The origin of the reference frame is taken at the center of the charged particle and r = r ^r. 
Then, for £1 > (or £1 < 0), n tends to be parallel (or perpendicular) to f near the charged 
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particle. For ei > the decrease of the electrostatic energy is estimated as 

Z2e2 / 1 1 



where the radius R plays the role of the lower cutoff length and the upper cutoff length 
£ is assumed to be longer than R and shorter than the screening length A. For £i < 0, 
£i in (21) should be replaced by |ei|/3 (because the angle average of (e^ • f)^ is 1/3). In 
the literature of physical chemistry, much attention has been paid to the decrease of the 
electrostatic energy around a charged particle due to polarization in the surrounding fluid 
(solvation free energy). Its first theoretical expression is the Born formula AE = Z^e^(l/eo — 
1)/2R}''^'^^ The additional orientation degrees of freedom in liquid crystals yield (21). 

We are assuming that n is appreciably distorted from Bz in the space region R <r <£. 
The Frank free energy is estimated as 

Fo ~ 7tK£. (22) 
We determine i by minimizing Fq + AF^ to obtain 

e = {\ei\/2T:elKf/'^Ze 

= {\e^\/2T,eof'\a£^fl^Z. (23) 

Here = /sqT is the Bjerrum length and we have set K = k^T/a with a being a micro- 
scopic length. The condition of strong orientation deformation is given by < ^ or 

R/Z <{\ei\/2T,eofl\atjif'\ (24) 

If this condition does not hold, the cfi"cct of the electric field on the ncmatic order becomes 
weak. Here £b ~ 500A for Cq ~ 1 ^-nd T 300K. Therefore, for a microscopic radius R, (24) 
can well hold unless |ei|/eo is very small. For salt with small concentration consisting of ions 
with charge Z^ and density pa (a = 1,2,...), the screening length is given by the Debyc-Hiickel 
expression = Aire^ Yla ^aPct/ ^O^^T . The criterion (24) is meaningful only for £ < A. If 
R < X < £, the director orientation is deformed in the space region R < r < X and l/£ in (21) 
should be replaced by 1/A. 

To illustrate the deformation of n around charged objects in equilibrium, we have numer- 
ically solved 



n X 



= 0, (25) 



and (10) in two dimensions by assuming n = (cos ^, sin ^) (or riz = 0). A charge is placed in 
the hard-core region (x^ + y^)^/^ < R. In three dimensions this is the case of an infinitely 
long charged wire with radius R and charge density a, in which all the quantities depend 
only on x and y. The solution can be characterized by the three normalized quantities, £i/£o, 
a* = a I {eqKYI'^ ^ and R/£-b- Here we set a* = 2.4 and R/£-b = 2. We discretize the space into 
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a 200 X 200 lattice in units of under the periodic boundary condition in the x direction, 
so the system width is L = 200^b- The electric potential vanishes at y = and L. In Fig. 2, 
n tends to be parallel to the y axis far from the origin and the spacing between the adjacent 
bars is 2£b- The director is in the radial direction for ei/eo = 0.4 in (a) and is perpendicular 
to r for ei/eo = —0.4 in (b) near the origin. A pair of defects are aligned in the x direction 
in (a) and in the y direction in (b). Similar defect formation was numerically realized in two 
dimensions for a neutral particle with the surface interaction (26) below. ^' 

5. Charged colloidal suspension 

The nematic orientation around a large particle without charges^' is determined by 
competition of the Prank free energy Fq and a microscopic anchoring interaction at the particle 
surface expressed as 

F^ = ~W^JdS{n-rf. (26) 

Here J dS is the surface integration, and r is the normal unit vector at the surface. The degree 
of anchoring is represented by the dimensionless parameter, 

Ma = W^R/K. (27) 

In charged colloidal suspensions, the distortion of n due to the surface charge can be more 
important than that due to the anchoring interaction (26). Note that (24) is well satisfied 
even for large R if the ratio Z/R can be increased with increasing R. Notice that the ionizable 
points on the surface is proportional to the surface area AtiR? . However, the problem is very 
complex, because the counterions themselves can induce large deformation of the nematic 
order (because of their small size) and tend to accumulate near the large particles. 

For simplicity let the deformation of n around the counterions be weak. Furthermore, we 
assume that the screening length A is shorter than R and the nematic order does not vary on 
the spacial scale of A except for the defect-core regions.® Then is approximated from (20) 
as Fe ^ -IFe / dS{n ■ rf /2 with 

We = eiZ'^e^X/A'KelR'^. (28) 

This is of the same form as Fa in (26) and the strength of anchoring is represented by /Xes = 
{Ws, + We)R/K. Here the Debye-Hiickel screening becomes extended in the dilute limit. In 
such cases and in the limit of large Z/ R, the electric field around a large particle decays on 
the spatial scale of the Gouy-Chapman length, 

A = k^T/eE^ = R^/i-BZ, (29) 

where = Ze/ 'ii^eoR^ is the electric field at the surface. If l^ZjR 2> 1, we surely have 
A < i?. Then we obtain We = (£ifeBr/47r£o)-^/-R^ and 

/xeff = (Wa/K)i? + {eikBT/4.'KeoK)Z/R. (30) 
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6. Concluding remarks 

In Section 2, we have presented a scheme of Ginzburg-Landau theory for the electric field 
effects in liquid crystals, where the dielectric tensor depends on the orientation order as in (3). 
Generalizations to more complex situations such as ferroelectric cases are straightforward. In 
Section 3, the director fluctuations have been examined in applied electric field for arbitrary 
strength of the dielectric anisotropy. It is desirable if (17) could be confirmed by light scattering 
experiments in systems with not small ei/eo- In Sections 4 and 5, the distortion of the nematic 
order around a charged particle has been studied for the first time. The condition of strong 
orientation deformation due to the charge effect is given by (24) if the charge is not screened 
in the range r < i with ^ being defined by (23). In colloidal suspensions the presence of the 
counterions makes the problem very complex, where we have obtained (30) for the case in 
which the colloid surface charges are screened within a thin layer with thickness \ <^ R. 

We finally mention a similar charge effect recently predicted for a near-critical polar binary 
mixture.^) In such fluids the dielectric constant eic) strongly depends on the composition c and, 
as a result, the phase separation behavior is strongly affected by doped ions. In particular, 
if charged colloid particles are suspended, a charge-density-wave phase should be realized 
for (|£i|/4£o)^/-R > (ttCo/^b)''^'^^ at low temperatures, where £i = d£{c)/dc and Co is a 
microscopic wave number (=the coefficient in the gradient free energy). Here large size of 
{\£i\I£q)Z/R is required as in (24). 

This work is supported by Grants in Aid for Scientific Research from the Ministry of 
Education, Science, Sports and Culture of Japan. 

Appendix 

We connect two capacitors in parallel as in Fig. lb. One contains an inhomogeneous dielectric 
material under investigation and the other is a large capacitor serving as a charge reservoir. 
The area and the charge Qq of the large capacitor are much larger than S and Q of the 
smaller capacitor, respectively. We are supposing an experiment in which the total charge 
Qtot = Qo + Q is fixed and the potential difference is commonly given by $ = Qo/C'o) where 
Co is the capacitance of the large capacitor. Obviously, in the limit Q/Qq ~ S/Sq 0, 
the deviation of $ from the upper bound $tot = Qtot/C'o becomes negligible. Because the 
electrostatic energy of the large capacitor is given by Eq = Q'^/2Cq = (Qtot 72^*0) (l — Q/Qtot)^; 
we obtain 

Eo = (Qtot/2C^o) - ^Q, (A.l) 

where the first term is constant and a term of order (Q/Qtot)^ is neglected. Therefore, for the 
total system including the two capacitors, the relevant free energy is given by G in (11). 
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Fig. 1. (a) System of a capacitor and an inhomogeneous dielectric material at fixed capacitor charge 
Q. The potential difference $ is a fluctuating quantity, (b) Two capacitors connected in parallel 
with charges Q and Qo- The smaller one contains an inhomogeneous dielectric material, and the 
larger one a homogeneous dielectric material. In the limit Q/Qo 0, the potential difference $ 
becomes fixed, while Q is a fluctuating quantity. 
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Fig. 2. The director in two dimensions for ei/eo — 0.4 in (a) and for ei/eo — —0.4 in (b) near the 
origin. It tends to be along the y axis (in the vertical direction) far from the origin. The spacing 
between the adjacent bars is 2£b- 
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